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Thermal fluctuations provide the main source of large scale density perturbations in warm infla- 
tionary models of the early universe. For the first time, general results are obtained for the power 
spectrum in the case when the friction coefficient in the inflaton equation of motion depends on 
temperature. A large increase in the amplitude of perturbations occurs when the friction coeffi- 
cient increases with temperature. This has to be taken into account when constructing models of 
^\ ' warm inflation. New results are also given for the thermal fluctuations in the weak regime of warm 

' inflation when the friction coefficient is relatively small. 
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I. INTRODUCTION 

T-H ■ Observations of the cosmic microwave background have provided strong evidence that the primordial density fluc- 
04 ] tuations have a power spectrum very close to that predicted by inflationary models of the very early universe 

According to the standard theory, these density fluctuations can arise from quantum vacuum fluctuations during 
inflation Another possibility, which also fits the data, is that thermal fluctuations provide a source of density 
fluctuations Q . This is what happens in the warm inflationary scenario 0] . 

During an inflationary era particles are being produced continually but their density is also rapidly diminished 
, by the expansion of the universe. What distinguishes warm inflationary models is that this particle production is 
Qh' sufficiently strong, compared to the effects of expansion, to produce a non-negligible particle density. Warm inflation 
Q I occurs when the quantum statistical fluctuations in the particle number are large enough to influence the classical 
^ inflaton field and produce density fluctuations [1, [§] . 

The description of coupled fluctuations between classical sources and quantum radiation flelds can be traced back to 
Ci . Einstein's classic discussion of a mirror immersed in black body radiation ■ Stochastic fluctuations in the radiation 
field lead to fluctuations of the mirror. In an analogous way, fluctuations in the radiation field are transfered into 
I inhomogeneities in the inflaton field. In a non-expanding universe, standard thermal field theory methods can be 
J> ■ used to reduce the description of the system to a pair of coupled stochastic equations, one for the inflaton and one for 
&^ [ the radiation density [ll|, |T2| . The simplest case to analyse is where the radiation is close the thermal equilibrium, 
. and then the particle production rate, the dissipation in the inflaton equation of motion and the source terms in the 

inflaton equation are all related by a single friction coefficient T [TT], [iM, [H, [3 • 
. There has been progress in extending the theory of thermal fluctuations to an expanding universe , but it seems 
VO ' that the basic equations which we require can be determined from the equivalence principle [1, [l^ . This approach 
, will be pursued further in Sect. Ill, where we present some new modifications to stochastic source terms in curved 
space. Our results, however, are broadly in line with those presented elsewhere which imply a close to scale-free power 
spectrum of density perturbations with small corrections depending on the slow-roll parameters [13, EM [l^ ■ 

The main objective in this paper is to present results for the power spectrum of fluctuations in inflationary models 
where the friction coefficient F depends on the temperature T. Temperature dependence seems to be a common 
^ feature in explicit calculations of the particle production rates. In particular, there are high-temperature models of 
d ' warm inflation with friction coefficient F (x [l^, [2lj and the two-stage decay models [22] with friction coefficient 
F oc [20|. High-temperature models are difficult to realise without producing large thermal corrections to the 
potential, which make the models unstable [23| . and most consistent warm inflationary scenarios use the two-stage 
decay mechanism (see Ref. Q for a review). 

Individual models with temperature dependent friction coefficients have been examined using the numerical code 
written for Ref. [s!], which solves coupled perturbation equations for the time evolution of the infiaton, radiation and 
metric. The model with F oc described in Ref. Q produced an oscillatory power spectrum, but it has not been 
possible previously to formulate general laws for the functional dependence of the power spectrum. 

An analysis of the coupled perturbation equations is presented in Sect. IV. We have been able to find approximate 
solutions to the perturbation equations on sub-horizon scales to determine general formulae for the power spectrum. 
These formulae have been checked against a numerical solution of the stochastic perturbation equations in section V. 
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The most important new results apply to the 'strong' regime of warm inflation where T ^ 3H and H is the 
expansion rate. There is a large enhancement of the density perturbation amplitude when T (x T^, and a large 
diminution when T cx T^^. In practice, the power spectrum of density fluctuations is fixed by observation, and so 
the enhancement of the amplitude is equivalent to a reduction of the inflationary energy scale, and visa versa. We 
have also found that the oscillations of the power spectrum found in Rcf. [8] are removed by the ensemble averaging. 
These results are an essential ingredient for the construction of warm inflationary models. 



II. WARM INFLATION 



We shall assume that particles are produced during the inflationary era and that the particle interactions are strong 
enough to produce a thermal gas of radiation with temperature T. The necessary conditions for this to happen where 
discussed in Ref. • Warm inflation is said to occur when T is larger than the energy scale set by the expansion rate 
H. The particle production is associated with a dissipative effect on the inflaton, whose equation of motion becomes 

in 

0+ (3i7 + r)0 + y^ = (1) 

where r(0, T) is a friction coefficient, H is the Hubble parameter and is the (f> derivative of the inflaton effective 
potential V^(0, T). The effectiveness of warm inflation can be parameterised by a parameter r, defined by 

When r ^ 1 the warm infiation is described as being in the strong regime and when r ^ 1 the warm inflation is in 
the weak regime. 

Consistent models of warm inflation [2^ require a suppression of thermal corrections to the inflaton potential, so 
that the effective potential separates into inflaton and radiation components 

y(0,T) = F(0) + p,(r), (3) 

where pr is the radiation density 

Pr = ^9*T\ (4) 
In this case, the time evolution is described by the equations 

^+i3H + T)<P + V^^0, (5) 

Pr+4.Hpr=T4,^, (6) 

SH^ = A-kG {2V + 2pr + (7) 

During inflation we apply a slow-roll approximation and drop the highest derivative terms in the equations of motion, 

m{l + r)^ + V^^Q, (8) 
AHpr = T^\ (9) 
3i?2 St:GV (10) 

The validity of the slow-roll approximation depends on the slow roll parameters defined in Q, 



The slow- roll approximation holds when e<Cl-l-?', Ty^l + r and /3 <C 1 + r. Any quantity of order e/(l -\-r) will be 
described as being first order in the slow roll approximation. 
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III. INFLATON FLUCTUATIONS 



Thermal fluctuations are the main source of density perturbations in warm inflation. Thermal noise is transferee! to 
the inflaton fleld mostly on small scales. As the comoying wavelength of a perturbation expands, the thermal effects 
decrease until the fluctuation amplitude freezes out @. This may occur when the wavelength of the fluctuation is 
still small in comparison with cosmological scales. 

In this section we shall review the equation for the inflaton fluctuations in the case where the friction coefflcient is 
independent of temperature, i.e. 

T^m. (12) 

We shall see how the stochastic inflaton dynamics transfers to an expanding universe and obtain new results for the 
case r < 3H. 

The behaviour of a scalar field interacting with radiation can be analysed using the Schwinger-Keldysh approach to 
non-equilibrium field theory [2^, [2^. In fiat spacetime, when the small-scale behaviour of the fields is averaged out, 
a simple picture emerges in which the field can be described by a stochastic system whose evolution is determined by 
a Langevin equation 27[. This takes the form 

- V^^ix, t) + T^{x, t) + V^ = (2rT)i/2^(a:, t), (13) 

where is the flat spacetime Laplacian and f is a stochastic source. For a weakly interacting radiation gas the 
probability distribution of the source term can be approximated by a localised gaussian distribution with correlation 
function [lll,[ll, 

{ax,mx',t'))^5^^\x-x')S{t-t'). (14) 

We shall restrict ourselves to this gaussian noise approximation. 

We can use the equivalence principle to adapt the flat spacetime Langevin equation to an expanding universe during 
a period of warm inflation when T > H. The Langevin equation will retain its local form as long as the microphysical 
and thermal scales in the problem are small compared to the cosmological scale 0, [31 ■ However, the rest frame of 
the fluid will have a non-zero 3— velocity with respect to the cosmological frame and we must include an advection 
term. The Langevin equation becomes [l6| 

4>{x, t) + (3i7 + r)0(x, t) + ra-2<5a0(x, t) + V^- a-~^d^(t){x, t) = {2T^fiTf'^£,{x, t) (15) 

where a is the scale factor and is the Laplacian in an expanding frame with coordinates x". The correlation 
function for the noise expressed in terms of the comoving cosmological coordinates becomes, 

{£,{x,t)^{x' ,t')) ^ a-'^(2Tif5'^'^\x - x')5{t ^ t'). (16) 

A new parameter Fofi has been introduced because at this stage we cannot determine the effects of the expansion on 
the noise term. 

Next we treat the source term as a small perturbation and expand the infiaton field 

(l){x,t) -> (l){t)+5(j){x,t) (17) 

where 5(j) is the linear response due to the source ^. This expansion is substituted into the Langevin equation and 
then we take the Fourier transform, 

(5<ii(k, t) + [iH + F)(5</)(k, t) + V00(5(/)(k, t) + fc2a-2(50(k, t) = (2Feffr)i/2^(k, t). (18) 

The infiaton will also generate metric inhomogeneities, but with a suitable choice of gauge, these can be discarded 
on sub- horizon scales (see [l^). We shall use a uniform expansion rate gauge. Eq. (fT8|) applies on scales which are 
intermediate between the thermal averaging scale and the horizon scale. Later, we use a matching argument to extend 
the fluctuations to large scales. 

The analysis of the Langevin equation can be simplified by introducing a new time coordinate z — k/{aH) and 
using the slow roll parameters (fTT|) . We are led to the equation for (5(/)(k, z), 

(1 - ~ef5(j)" - (3r + 2)(1 - ~e)z-^c^' + ~e'5<i)' + 3r/z-2j0 + 5</) ^ {2T,sTY'^{1 ~ ef'^S, (19) 
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where a prime denotes a derivative with respect to z and e = e/(l + r). This reduces to 

(50" - (3r + 2)z-M0' + ^0 = (2roffr)i/2^~ (20) 

when we keep only the leading terms in the slow roll approximation. The noise term has been rescaled so that the 
correlation function is now 

(C(k, z)e>', z')) = (2^)3^(3) (k - k')(5(z - (21) 

The perturbation equations can be solved using Green function techniques. The solution is 

> 

G(z, z'){z'f-''- (2reffT)i/2^(z') dz'. (22) 

The retarded Green function G(z, z') can be found in terms of Bessel functions, 

G(^z,z')^'^z''z"' {J,Az)Y^{z')- J,{z')Y,{z)). for z < z'. (23) 

where 

3, , r + 3iJ 
^=2(1 + ^) = ^^. (24) 

Corrections due to the time dependence of v are similar in size to terms which we have already discarded in the slow 
roll approximation. 

A. Inflaton power spectrum 

The inflaton power spectrum P^{k, z) is defined by 

((50(k, z)5<^(k', z)) = F^(fc, z)(2^)35(k + k') (25) 
Substituting the first order inflaton perturbation from (1221) and the correlation function (j2ip gives 

P^{k, z) = (2reffT) / dz'G{z, z'fiz'f-^'' (26) 

J z 

in the slow-roll approximation. 

Integrals of this type are examined in appendix [A] The large z approximation corresponds to length scales much 
smaller than the horizon size, and we can use this to examine the flat space limit, 

P^ik, z) ^ (2reffr) k-^^^^^ (27) 

If the equivalence principal applies, then the result must agree with a calculation using flat space thermal field theory 
(where the physical momentum is fca^^), 

P^(fc,z) = fca^^Tfc^^ (28) 

This allows us to read off the value of the constant Fcff from a comparison of Eq. (|?7)) and Eq. (^5]) using Eq. ([M]) , 

Teff = r + H. (29) 

The value of the coefficient is required for any analysis of the density perturbations in the weak regime(r < 1) of 
warm inflation. Such calculations have previously used the friction term F + 3-ff for Fcg, but fortunately this has 
little effect on the final power spectrum [3]. 

When r is small and u ~ 3/2, the perturbation amplitude at horizon crossing (z = 1) can be found numerically, 

P^{k, I) M 1.526 k-'^ HT (30) 
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When r is large, the perturbation amphtude at horizon crossing is governed by a saddle point in the integral (j26p . 
given according to eq. (|A11[) by, 

k^a^iHitpntF))'/^ (31) 

We call the time in eq. (|3ip the freezeout time tp for the mode k. The freezeout time always precedes the horizon 
crossing time. An analytic approximation to the integral valid for large values of r is given in eq. (|A5p . With this we 
recover the result derived in 8] , improving on an earlier results in Q , 

P4k,l)^k-^^{THy/^T. (32) 
An approximation which works remarkably well for both large and small values of r is given by 0] 

P^(fc,l)«A:-3^^i/i/2(r + 3i/)i/2r. (33) 
We shall compare the accuracy of this approximation to numerical results in Sect. V. 

IV. COUPLED FLUCTUATIONS 

In this section we analyse the density perturbations in the case where the friction coefficient in the inflaton equation 
depends on temperature, i.e. 

r = r(T,0). (34) 

Temperature dependence in the friction coefficient is a feature of many models [20l |. and numerical calculations have 
shown that this can lead to interesting effects on the density perturbations (sj. The temperature dependence of the 
friction coefficient is parameterised by a coefficient c, 

c=^^ (35) 

The consistency of warm inflation requires only that c < 4, and values of c = 3 and c — ~1 arise in calculations of 
the reheating term for different models. 

The non-zero value of c leads to additional terms in the fluctuation equations through 

sr = TtSt + r^S(j). (36) 

The extra S(j) terms depend on the slow-roll parameter f3 and they are first order in the slow-roll approximation, but 
the 6T terms are leading order and couple the perturbations of the inflaton and the radiation. 

A. Radiation fluctuations 

There are two sources of fluctuations in the radiation. The first source is purely statistical and caused by the 
microscopic particle motions. These fluctuations can be described as photon number fluctuations Q , thermodynamic 
or quantum field fiuctuations Q, but all approaches agree on the power spectrum Ps for Sp,. in the flat space limit, 

P, = AprTa-^. (37) 

An estimate of the relative size of these fluctuations can be made by comparing the statistical fluctuations to the 
scalar field potential fiuctuations V0(50 calculated using Eq. ([55)1 at horizon crossing (fc = aH), 



vlP^ y2^i/2(r + 3i/)i/2rfc-3 (i + r)3/2- 



These stochastic radiation fluctuations are relatively unimportant in both the weak (r ^ 1) and strong (r ^ 1) warm 
inflationary regimes. 
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The second source of fluctuations, which turns out to play the dominant role, is due to the inhomogeneous energy 
flux from the inflaton field. We can describe these fluctuations using first order perturbation theory. This can be 
formalated in many equivalent ways [1, [11] , but the formulation we use here is based on Ref. [2^ , and Ref. [l^ . The 
transfer of momentum and energy into the radiation is described by a momentum flux J and an energy flux JQ, given 

by 

J = ~T(i)S(j), (39) 
6Q ^ ^"^ST (40) 

If we keep only the leading terms in the slow roll approximation, the metric perturbations drop out and the energy 
density perturbations in a uniform expansion gauge satisfy the equation 

t2 , 1 ,2-2 \ x„ _ u2-2 



Spr + 9H6pr + \^20H-' + ^k^a-^j Spr = fc^a" V + 5HSQ + 5Q (41) 

The energy flux (j40p can be expressed in terms of Spr using Eqs. (|35p and p6p by 

SQ = cHSpr, (42) 

The radiation fluctuations therefore satisfy 

6pr + (9 - c}HSpr + (^{20 - 5c)H^ + ^k^a'^^ 5pr = -k^a-^{T^) Scj) (43) 

The equation for the inflaton fluctuations can simply be obtained by adding a term dT(f> to Eq. (jl8p , 

64> + m{l + r)S<j3 + k^a-^S(j3 + Scr(r<j3)-^Spr = {2r^sTy^^^ (44) 

There are now two coupled equations which we can use to determine the fluctuation power spectrum P^. 

An important feature of the coupled system is the growing mode d4> oc a^'^ present for times such that k ^ aH ^ k/r. 
This replaces the constant mode responsible for freezing out the perturbation amplitude on sub-horizon scales for the 
uncoupled system. The concept of 'freezout' is a useful one nevertheless, as it defines the time when the noise terms 
leave their imprint on the final fluctuation amplitude. We will retain the terminology, but with the understanding 
that 'freezout' now means approaching the growing mode. 



B. Asymptotic expansions 



We shall determine the fluctuation amplitude for large r by using a matched asymptotic expansion technique. To 
begin with, we convert the fluctuation equations to the variable z — k/{aH) and drop the terms which are flrst order 
in the slow-roll approximation, 

6(f>" - {3r + 2)z~^6(j)' + 6(f> + crz-^{r(j)y^Spr = {2r^tfT)^/'^i (45) 
dp';- {8 -c)z-^Sp'^ + {20 -5c)z'^6pr + \spr + {r(j))S(P ^ (46) 

We suppose that r is a large parameter and cr — 0(r). 



1. The range z « r^^"^ 

In the range z « r^l"^ ^ the stochastic source terms are relatively unimportant and the leading order contribution 
to Eq. (|45|) for large r is simply 

5pr « c-\T^)z—6(j). (47) 
oz 



Substitute this into Eq. to get 
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The solutions to this equation are Meijer-G functions, (see |30[), 

r2 



z' 
12 



1 - 3c/2 
2-c/2 5/2 



G 



30 ( „±i7r • 
13 



12 



1 - 3c/2 

2-c/2 5/2 



(49) 



We are interested in the solution which grows most rapidly with time. Since the time evolution proceeds from large 
to small z, this is the solution which decays most rapidly with z. 

The asymptotic expansion of the Meijer-G functions for large z can be found using standard techniques, in particular 



Vl2 



1 - 3c/2 

2-C/2 5/2 



-3c/2 



rflX3c/2)r;^.(3c/2 + 5/2)r;j,(2 + c) 1 



9c(3c + 5)(c + 2) 



where V ^ is the Gamma function. When c is positive, this is growing mode cx a'^'^. For the other solutions 



G 



30 ( „±i7r ■ 
13 



12 



1 - 3c/2 

2-C/2 5/2 



(3+c)/2 



,=Fjz/\/3T«'r(c-l)/2 



These oscillatory solutions grow more rapidly with time than the power-law solutions when c < —3/4. We set 



AG 



1 - 3c/2 



TIBG 



2~cl2 5/2 
where A and B are constants and TZ indicates taking the real part 



'30 ( „i7r ' 
13 



12 



1 -3c/2 

2-c/2 5/2 



(50) 



(51) 



(52) 



2. The range 1 << z « r 

In the range 1 << z << r we have to take into account the noise term. If c > —3/4, then Spr has power-law 
behaviour and the leading terms in Eq. (|46p give, 

6pr ~ -3(r0)(5</). (53) 

Substitute into Eq. P5|) . 

S(t>" - (3r + 2)z-i(5(/)' + S(t>- 9crz-^6<j) = {2TTy/^^ (54) 
As before, we solve this using a Green function, 

G'^iz, z) = 'i^z^z'" {Ju+3.c{z)Y^+y,c{z') - Jy+3c{z')Y^+y.c{z)) for z < z . (55) 

where v = 3(1 -I- r)/2. The power spectrum is given by 

dz' G^{z,z'f z'^-^P (56) 

If we further restrict to 1 << z << r^/^, then the integral can be approximated (see appendix A). Comparing to the 
previous result for c = 0, we get 

Pf°(fc,z) ^ rH(3c + 3/2) /3y- 

P,-0(fc,l)^ VR{i/2) {z^J • ^"'^ 

This result is not valid at horizon crossing (z — 1), and so a matching procedure is required. When c < —3/4 the 
solutions oscillate and Eq. ([55)1 is no longer valid. The best we can do with a simple analysis is to conclude from Eq. 
dSI]) that Pf°{k,z) oc (z^/r)'^+'^, multiplied by the ensemble average of oscillatory terms. We shall investigate this 
case numerically in the next section. 
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3. Matching in the range 1 << z « r^^"^ 

For c > —3/4, both approximations are valid (and agree) in the the range 1 << z << r^l"^ . Comparing 84? from 
Eq. dsn) and from gives B = and 



We are able to use Eq. ((52)) at z = 1 



P;^°(fc,l) I ( 1^ 



1 - 3c/2 

2-C/2 5/2 ; [ rfl(3c + 3/2)^3. 



P^^=0(fc,l) 23-rfl,(3c/2)rfl(3c/2 + 5/2)ri?(2 + c) rfl(3/2) 
Numerical values of the result for different values of c lie close to a simple power-law approximation 



(59) 



, 3c 

r 



(60) 



where Tc varies slowly in the range < c < 3, for example r\ k, 8.53, ri k, 7.66 and r-^ k, 1 .11 . 

C. Damping effects 

Damping effects can suppress the radiation fluctuations and may be important when the thermal equilibrium of the 
radiation is imperfect. In warm inflation thermal equilibrium is disrupted by particle production from the evolving 
inflaton field which competes with the thermalising effects of self-interactions within the radiation gas The 
effectiveness of the thermalisation process can be parameterised by the value of iJr^, where is the relaxation time 
of the radiation. 

The most important damping effect in a pure radiation gas comes from the shear viscosity r\ (3l| . which modifies 
the perturbation equation (|4ip to, 

8pr + (9iJ + r]k^a-^)Spr + (lOH^ + ^fc^a^^ ] Spr ^ k'^a-^J + 5HSQ + SQ, (61) 



where [311 



3 



V = T7Pr-Tr. (62) 

15 



The power-law behaviour of the growing mode is affected by the new damping term, so that a similar analysis to the 
undamped case now gives 



-1^, / ^ \ 3c' 



prik,i ) 

p^^ik, 1 

where 

2c 



l + (l + 24ci/Tr/5)i/2' 



(63) 



(64) 



The equilibrium is maintained when Hrr <C 1 and in this case c' ~ c. However, the value of c' is still significant even 
when HTr ^ 1 and the thermal equilibrium is very marginal. Enhancement of the fluctuations therefore appears to 
be a fairly robust phenomenon. 
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V. NUMERICAL SIMULATION 



The full set of perturbation equations for warm inflation, including the metric perturbations, where solved numeri- 
cally by Hall et al. in ref. 0]. In the case where the parameter c = — 1, the power-spectrum of the perturbations had 
a series of oscillations superimposed on the usual power-law behaviour in k. 

The computer code used in ref. [1] uses the stochastic equations to set up the initial conditions and then evolves 
the full perturbation equations deterministically. The drawback with this approach is that the code only runs once 
for each mode, and does not necessarily produce the same results as an ensemble average of the stochastic system. 
In this section we shall solve stochastic equations for the perturbations, but we work to leading order in the slow-roll 
approximation and neglect the metric perturbations. The deterministic code helps provide solid evidence that the 
metric perturbations are small on sub- horizon scales, confirming the same conclusion obtained analytically in ref. p^ . 

We shall use rescaled variables 

y = fc3/2(2r,gT)~i/2j0, w = 0^{2T,sT)-'/^{T^r'Spr. (65) 
The perturbation equations Eqs. (l45l) and (|46p become 

y" - (3r + 2)z-'^y' + y + crz-'^w = | (66) 
w" ~ (8 - c)z-'^w' + (20 - 5c)z-^w + \w + y = (67) 

where ^ is a normalised Gaussian random variable. The scalar field power spectrum is given by 

- (2reffr) k-^ {y{zf) (68) 

Departures from a scale-free spectrum arise from the dependence of F, H and T on the fc-dependent horizon crossing 
(or freezeout) time. This can be determined from a numerical solution of the background homogeneous field equations. 
Our aim here is to examine how {y{z)'^) depends on r and c. 

Numerical results for c — 0, c — 1 and c = — 1 are shown in figures [TJH These have been obtained using a fourth 
order Runge-Kutta integration scheme and the Box-MuUer algorithm for generating Gaussian random variables. Fig. 
[Ushows the time evolution of the inflaton field, plotted against z = —k/{aH) for convenience. The histogram on the 
right shows the distribution of field values at horizon crossing z = — 1 and fits a Gaussian distribution. The other 
graphs show (y^) at horizon crossing for different values of r = T/3H. 

The curves superimposed on the data are based on the theoretical results obtained in the previous section. For 
c > —3/4, Eq. (pO]) when converted into the rescaled variables corresponds to 

{y')^\{iy^'r;''r''-'^'- (69) 
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The theoretical prediction is excellent for c = and r ^ 1. The data can also be approximated by rational fractions 
over the entire range, for example the approximation 



(y^) - ^ii±l)''' (70) 

works extremely well for c = and leads to the prediction for the inflaton power spectrum given in Eq. (|33p . 
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FIG. 2: The variance of the inflaton at horizon exit from the numerical simulations has been plotted against r = T/3H for 
c — {r independent of T). The data are averages over 1000 runs. The curve shows the approximation Eq. (|70|l . 



The data for the c = 1 case have been fit to a curve of the form Air^^^ + Bir^^'^ based on the theoretical prediction. 
The need for next-to- leading order corrections to the asymptotic formula is indicated by the analysis, for example in 
Eq. ([50)1 . The coefficients for the best fit are shown in table HI and correspond to = 8.2 in Eq. ([5^ . which is 5% 
from the theoretical prediction — 8.5. Results for c = 2 and c = 3 using the function Acr^'^~^/^ + Bcr^'^~^/'^ are 
also shown in the table. 

The most significant feature in the data for the c = — 1 case is that we have found no evidence for the oscillations 
in the amplitude found previously in Ref. @ . This is presumably due to the fact that the numerical code used in Ref. 
Q only performs one run per mode, whereas we average over many runs. The results have been fit to a curve of the 
form A-ir~^ + B_ir~^, in accordance with the leading order (z^/r)'^+'^ behaviour of Scfy^ suggested in the previous 
section. The coefficients for the best fit are given in table HI 



TABLE I: CoefBcients in the numerical fits for the variance of the inflaton at horizon exit with their range of validity. 



c Ac Be range 

3 4.51 X 10"^ 3.34 x lO"'^ r > 100 

2 1.22 X lO"'' 4.39 x 10"" r > 20 

1 4.66 x 10"" 2.94 x 10"^ r > 10 

-1 5.30 -37.1 r > 20 
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FIG. 3: The variance of the inflaton at horizon exit from the numerical simulations has been plotted against r = F/Si? for 
c = 1 (F oc T). The data are averages over 1000 runs. The curve shows the numerical fit A\r^^'^ + B\r^^^ . 




FIG. 4: The variance of the inflaton at horizon exit from the numerical simulations has been plotted against r = T/3H for 
c — —1 (F cx T^^). The data are averages 1000 runs. The curve shows the numerical flt A-ir~^ + B-ir~^. 



VI. DENSITY FLUCTUATIONS 



We have seen that the small scale inflaton fluctuations freeze out by the time they cross the horizon. Whilst these 
fluctuations are freezing out, the metric fluctuations are relatively small (in the uniform expansion-rate gauge, for 
example). Eventually, the wavelength of the perturbations becomes far larger than the effective cosmological horizon 
and the metric perturbations become important. On these large scales it becomes possible to use a small-spatial- 
gradient expansion (first formalised by Salopeck and Bond [l^l). This approach allows us to define the curvature 
perturbation C so that it is conserved even in the non-linear theory [33, 34, 35]. The spectrum can be obtained, to 
high accuracy, by matching the small and large scale approximations at horizon crossing. 

The fluctuations can be described entirely by the conserved expansion fluctuation C on constant density hypersur- 
faces [3^, which is defined for general hypersurfaces by 



c4ln(l + 2^) + i/- 



dp 



where ip is the spatial curvature perturbation. After using the slow roll equations. 



1 



C= - ln(l + 2(^) + 



H 



(71) 



(72) 
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where H = H{(j)) and (p = (j>{(t)) are given from eqs. ((8l ll0|) . 

Consider a uniform curvature gauge ip = and C = C(0)- When the inflaton perturbations are expanded as before 
in eq. (fTT]) . we have 

C = C<pS^ (73) 

where (f) subscripts denote derivatives with respect to (j) and = H/(j). Hence the power spectrum of density 
perturbations is 

Pdk)=eMk) (74) 

In the previous section we evaluated the small scale inflaton perturbations in a uniform expansion-rate gauge, but 
on sub-horizon scales the metric perturbations are relatively small and the inflaton perturbations in the uniform 
expansion-rate gauge and the uniform curvature gauge are approximately equal. 

In the strong regime of warm inflation we can use Eq (|60p for the inflaton power spectrum when r > Vc and 
c > —3/4, to give 

This can be inverted, using the slow-roll equations (f8l fT0|) . to express the basic energy scales of the inflationary model 
in terms of the fluctuation amplitude, r and slow-roll parameters. The value of the potential when the mode with 
wave number k crosses the horizon is given by 

= ^1/12^1/4^-3/4 1^^^ ^ [^J^) (3-73 X IQi^GeV), (76) 

for r > Tc and c > —3/4. The c dependence depresses the energy scale of inflation relative to models with c = 0. 
In most models this leads to an increase in T/SH and T/H, making warm inflation stronger than in the c — case. 
When c < the energy scale of inflation increases and the warm inflation becomes weaker. On the other hand, if 
primordial tensor perturbations are seen in the cosmic microwave background, this would favour a large energy scale 
for inflation and c < 0. 



VII. CONCLUSION 



We have analysed the thermal fluctuations present in warm inflationary scenarios to obtain results for the primordial 
power spectrum using both analytical techniques and numerical simulations. For the first time it has been possible 
to present general formulae for the power spectrum in models where the fiction coefficient in the inflaton equation of 
motion depends on temperature. These results are crucial for building models of warm inflation which are consistent 
with cosmic microwave background observations. 

The results for the primordial power spectrum of curvature fluctuations are as follows. For a friction coefficient F 
which is independent of temperature, the familiar approximation 

= ,77, 

fits both analytic and numerical calculations remarkably well for all values of F. This result can be applied whenever 
thermal fluctuations are the dominant source of curvature fluctuations, which is typically when T > H . 

The power spectrum is given by a different formula when the friction coefficient depends on temperature. For 
F cx T^, and r — T/3H ^ 1, the amplitude is enhanced when c > and reduced when c < 0. For accurate results we 
recommend using 

P^ik)^6^f{r)k~^ (78) 

where the function /(r) ~ AcV'^'^'^^^'^ + iJ^r'^'^"^/^ for c > and the coefficients are given in table HI 
For less accurate work, a reasonable approximation for c > is provided by 

Pdk) = ^^(1 + ry/' (l + f ) (79) 



13 



with ra « 7.3. This formula is consistent with the old result given by Eq. (|77p for c = and approaches the numerical 
results for large r. 

The most promising models of warm inflation to date use the 'two-stage' decay mechanism, for which the friction 
coefficient has c = 3 [33,[3l|. The new power spectrum is very different from the old result Eq. ([77]) . and therefore these 
models must be analysed using the new power spectrum. The effect of the new formula is to reduce the inflationary 
scale, and it is likely that the restrictions on warm inflationary models will be relaxed. 

It is known that warm inflationary models produce a significant amount of non-gaussianity, possibly enough to be 
observed in microwave background experiments [16, 39]. So far, the analysis of non-gaussianity has only been done 
in cases where the friction coefficient is independent of temperature. However, it seems possible that the analysis 
of the coupled perturbation equations initiated here could be extended to second order to obtain results for the 
non-gaussianity when the friction coefficient is temperature dependent. 



APPENDIX A: INTEGRALS 



We begin with an approximation to the integral 

F{z)= dz'G{z,z')^{z''f-*'' (Al) 



where — 3{1 + r)/2 and the retarded Green function G is given in eq. (|23p . The leading terms for large ly and fixed 
z, come from 

F{z) « / J.iz'fz'^-^^dz' (A2) 

4 Jo 

This is a standard integral, 

where Tn is the gamma function. We also have 



Hence, 



A similar integral arises when we replace G with the Green function G"^ in ((55)l . and in this the case the integral, 
which we now call i^^(z), has an asymptotic expansion with leading term 



tr{c + 3/2) f2i,Y' rir 



in the range z <C v^^"^. 

A slightly more challenging problem is the integral 



/oo 
G{z,zrz"-"'9iz')dz' (A7) 

where g{z) is a smooth function. We proceed as above to get 



Fiz,g)^ —z^'^Y^izY J,{z'Yz"-"'g{z')dz'. (A8) 

4 Jo 

For large v there is a Debye approximation for the Bessel functions which is valid in the range Q < z < 

J^{z) - (27rj^tanha)^l/2g^(tanha-a) ^^g) 
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where vsecha = z'. The relevant part of the integral becomes 

I'OO 

F(z,5) « ^z2''r,(z)2^.2-2^' / (^osh a)^'^-'^ c^^^**"^""") ^(z') rfa (AlO) 

where vsecha = z'. We find that there is a saddle point in this range at the value of a corresponding to z = zp, 
where 

Expanding about the saddle point gives 



J^(^,5)~^32^5(rF) (A12) 

which agrees with our earlier result when g{z) = 1. This saddle point is responsible for the phenomenon of 'freezing 
out' of the thermal fluctuations. The value of z decreases with time and the fluctuations always freeze out before they 
cross the horizon at ^ = 1. 
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